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Thermodynamic analysis of parallel-serial reaction in complex oxide systems is carried out. The use of the 
Karapet’yants method to supplement the refined Hooke-Jeeves method has made it possible to determine the 
optimum range of compositions studied and the most rational temperature interval for heat treatment of coat¬ 
ings. This makes it possible to substantially decrease subsequent experimental search and reduce it to a verifi¬ 
cation of the optimum interval that is thermodynamically identified. 


The possibilities of chemical thermodynamics go beyond 
establishing the energy possibility of a specific reaction. In 
particular, a method for calculating an equilibrium phase 
composition in batch mixtures heat-treated at temperature 
not less than 1473 K is rarely used. This method makes it 
possible to predict the qualitative composition of high-tem- 
perature reaction products and estimate their quantitative ra¬ 
tios in target material, which is significant for research of 
crystallizing silicate glasses and substantially speeds up sub¬ 
sequent experimental search; in fact, this search becomes re¬ 
duced to a verification of the optimum range identified for 
initial compositions 

A calculation of this kind has been implemented with re¬ 
spect to non-fritted high-temperature glazes, whose structure 
and properties are formed in the course of melting of an ini¬ 
tial batch mixture and become fixed in the cooling of heat- 
treated articles. The calculation was carried our for composi¬ 
tions with fundamental differences in the content of the man 
components (Table 1). The mechanism of phase interactions 
and the properties of coatings based on them will be different 
as well; therefore, prediction of the equilibrium composition 
of reaction products is interesting precisely for such oxide 
systems. 

It should be noted that real batch mixtures contain not 
only individual oxides but also complex natural and chemi¬ 
cal compounds (pegmatite, chalk, barium carbonate, waste 
potassium-iron catalyst, etc.). Identification of their precise 
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thermodynamic characteristics is problematic; moreover, for 
such components as the potassium-iron catalyst it is totally 
impossible. Determination of these parameters based on in¬ 
crements [1] will produce a too big calculation error. In view 
of this, and also considering that calculations are to be ful¬ 
filled for high temperatures (900 - 1700 K), we believed it 
more reliable to analyze possible reactions between oxides 
making up initial materials, which is reflected in the list of 
reactions proposed. This solution appears even more sub¬ 
stantiated considering that some elements of this approach 
had been validated before [2]. 

The qualitative and quantitative oxide compositions of 
non-fritted glazes selected for the study make it possible to 


TABLE 1 


Oxide 



Weight content, %, in glaze 


3„ 

3 iv 

3 

J VI 

7 

^21 

Z 33 

Z 36 

7 

^37 

Si0 2 

57.8 

53.3 

56.4 

49.3 

50.8 

48.0 

51.2 

A1 2 0 3 

7.6 

9.7 

7.9 

20.0 

20.0 

18.6 

20.0 

CaO 

16.1 

13.2 

12.2 

3.0 

3.5 

6.0 

4.5 

MgO 

2.8 

2.8 

1.6 

0.3 

0.5 

2.2 

0.4 

BaO 

6.2 

6.3 

6.2 

5.5 

5.0 

6.5 

5.0 

ZnO 

- 

- 

- 

16.0 

16.0 

14.9 

14.5 

Na 2 0 

2.1 

2.0 

2.1 

1.9 

1.9 

1.9 

1.9 

K,0 

2.2 

2.0 

2.2 

1.0 

1.0 

1.0 

1.0 

Fe 2 0 3 

5.2 

5.2 

5.1 

0.3 

0.3 

0.3 

0.3 

Cr 2 0 3 

- 

- 

- 

2.7 

1.0 

0.6 

1.2 

Ti0 2 

- 

5.5 

6.3 

- 

- 

- 

- 
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Fig. 1. Equilibrium compositions of products of high-temperature 
reactions for mixtures of series 3 at 1500 K: 1 ) diopside; 2 ) hele- 
nite; 3 ) forsterite, 4 ) barium ferrite, 5) mullite. 


suggest that the most significant processes in batch mixtures 
will proceed according to the following reactions: 

2CaO + a-Al 2 0 3 + p-Si0 2 = 2CaO ■ A1 2 0 3 ■ Si0 2 ; (1) 


3a-Al 2 0 3 + 2P-Si0 2 = 3A1 2 0 3 ■ 2Si0 2 ; (2) 

CaO + MgO + 2Si0 2 = CaO • MgO ■ 2Si0 2 ; (3) 

MgO + Fe 2 0 3 = MgO • Fe 2 0 3 ; (4) 

2FeO + P-Si0 2 = 2FeO • Si0 2 ; (5) 

MgO + a-Al 2 0 3 = MgO ■ A1 2 0 3 ; (6) 

FeO + a-Al 2 0 3 = FeO • A1 2 0 3 ; (7) 

2MgO + P-Si0 2 = 2MgO • Si0 2 ; (8) 

ZnO + a-Al 2 0 3 = ZnO • A1 2 0 3 ; (9) 

2ZnO + P-Si0 2 = 2ZnO ■ Si0 2 ; (10) 

BaO + Fe 2 0 3 = BaO • Fe 2 0 3 . (11) 


TABLE 2 


Reac- 

Gibbs energy, kJ, at reaction 

temperature, 

K 

tion 

900 

1100 

1300 

1500 

1700 

1 

- 124.817 

- 127.571 

- 130.195 

- 132.748 

- 135.257 

2 

- 82.778 

- 82.867 

-82.881 

- 82.843 

-82.761 

3 

- 10.880 

- 17.330 

-56.867 

-68.567 

- 86.474 

4 

- 17.693 

-9.724 

+ 1.485 

+ 15.973 

+ 33.770 

5 

- 11.899 

-5.913 

-0.140 

+ 5.333 

+ 10.436 

6 

-52.135 

-61.206 

- 70.675 

- 80.454 

- 90.503 

7 

+ 3.954 

+ 17.431 

+ 29.785 

+ 40.654 

+ 49.748 

8 

- 57.803 

- 56.622 

-55.706 

-55.066 

- 54.694 

9 

- 14.429 

-98.167 

-205.154 

-333.833 

-483.085 

10 

- 18.652 

-9.701 

+ 1.207 

+ 13.893 

+ 28.235 

11 

+ 10.897 

-6.389 

-22.481 

-37.201 

-50.437 


G. V. Lisachuk et al. 

30- 


25- 2 



Z 21 Z 33 Z 36 Z 37 
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Fig. 2. Equilibrium compositions of products of high-temperature 
reactions for mixtures of series Z at 1400 K: 1 ) helenite; 2 ) mullite; 
3 ) gahnite; 4 ) sillimanite; 5) magnesium spinel; 6 ) forsterite. 

It should be noted that the oxidation of FeO to Fe,0 3 , 
2FeO + 0.5O 2 = Fe 2 0 3 , 

has been left outside the framework of the reactions consi¬ 
dered since it is hard to uniquely identify the oxide form of 
iron in the initial materials. Furthermore, the thermodyna¬ 
mics of this reaction does not have a fundamental effect on 
the equilibrium composition of the products of phase forma¬ 
tion due to subsequent formation of solid solutions with par¬ 
ticipation of both FeO and Fe,0 3 . Titanium and chromium 
oxide were not taken into account in the calculation either, 
since it is known that the former in an amount below 7% 
(here and elsewhere weight content) predominantly passes 
into the vitreous phase and the second one due to its high 
melting point will participate in reactions occurring at tem¬ 
peratures above the considered ones. 

The method proposed by Karapet’yants [1] was selected 
for solving this problem. Flowever, this method is not in¬ 
tended for solving a system of nonlinear high-order equa¬ 
tions expected in this problem, therefore, the Flooke - Jeeves 
method [3] was used in implementation with some necessary 
refinements [2], 

The temperature dependence of the Gibbs energy was 
preliminary determined for each reaction using the precise 
integration method [4], The results obtained with regard to 
the above considerations are listed in Table 2. 

Analyzing the obtained Gibbs energy values, it should be 
noted that the reaction of formation of hercynite (7) is ther¬ 
modynamically impossible within the entire temperature in¬ 
terval considered. Furthermore, three other reactions, namely 
(4), (5), and (10), have a restricted temperature interval of 
existence. 

The HUK-Ceramic program was used to calculate an ex¬ 
pected equilibrium composition of the products of the above 
reactions for each composition with a temperature variation 
step of 100 K. Based on the calculation data we found that 
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Fig. 3. X-ray diffraction pattern of coating 3 VI fires at 1553 K: 
O) barium ferrite BaFe 0 , 4 Fe 076 O 28g ; 0) barium ferrite BaFe 2 0 4 ; 
□) sanidine KAlSi 3 0 8 ; ■) diopside CaMg(Si0 3 ) 2 . 

the main products in series 3 is helenite 2CaO • A1 ? 0 3 • Si0 2 , 
whose content within the whole temperature interval is 
within limits of 20 - 29%. Diopside CaO ■ MgO ■ 2Si0 2 , for- 
sterite 2MgO ■ Si0 2 , and barium ferrite BaO • Fe 2 0 3 are also 
formed in perceptible quantities. 

The main reaction product in the composition series Z is 
3A1 2 0 3 • 2Si0 2 in an amount up to 25%; its accompanying 
components are helenite up to 15%, gahnite ZnO • A1 2 0 3 , 
and sillimanite A1 2 0 3 • Si0 2 . 

Figures 1 and 2 show comparative content diagrams of 
the main crystalline phases in the compositions considered 
depending on the formation temperature. 

The reactions of formation of fayalite 2FeO ■ Si0 2 , for- 
sterite 2MgO • Si0 2 , and decomposition of magnoferrite 
MgO • FcyO, occur intensely within a temperature interval of 
1000 - 1200 K in batch mixtures 3 0 and 3 IV . An opposite sit¬ 
uation is registered within the specified interval for composi¬ 
tions 3 VI , in which the content of magnoferrite increases in¬ 
significantly, and fayalite presumably decomposes, facilitat¬ 
ing the formation of barium ferrite. 

Phase transformations of fayalite and magnoferrite in 
all three mixtures end at temperatures of 1200- 1400 K, 
whereas the formation of barium ferrite continues and inten¬ 
sifies. Mixtures 3 0 and 3 IV are characterized by reactions of 
diopside with helenite, while the other phases have virtually 
constant contents. Iron oxide released after the completion of 
the reactions of formation of fayalite and magnoferrite is 
consumed to increase the content of BaO • Fe 2 0 3 and the for¬ 
mation of solid solutions of the type of (Ba, Fe)Al 2 0 4 or 
(Ba, Al)Fe 2 0 4 . The formation of color centers responsible 
for tinting of non-fritted glazes is possible on the basis of the 
compounds specified above. The series of phase transforma¬ 
tions stops within a temperature interval of 1400- 1600 K 
and the equilibrium compositions do not change. 

Analysis of estimated data for compositions of series 3 
suggests that coatings produced at a temperature of 1500 K 
will have different properties due to the presence of phases 
having an opposite effect on thermal stability (mullite and 



Fig. 4. X-ray diffraction pattern of coating Z 71 fired at 1523 K: 
■ ) gahnite; □) a-quartz. 

forsterite), chemical resistance (diopside and helenite), and 
color of coatings (iron-bearing solutions). 

The calculation of equilibrium compositions for series Z 
indicates that the formation of willemite in all of them stops 
at a temperature around 1400 K, which has a substantial ef¬ 
fect on the final composition of the reaction products, which 
exhibits a great variety. 

The temperature interval of 1200- 1500 K exhibits an 
intense reaction of the formation of zinc spinel. Such abrupt 
increase in its content in the mixture becomes possible due to 
the release of ZnO after decomposition of willemite and in¬ 
corporation of ZnO and Si0 2 into the glass matrix. 

Coatings of series Z will resemble each other in their 
main properties, since they contain a virtually constant quan¬ 
tity of gahnite, which is a compounds imparting high opacity 
(n = 1.82) and mechanical strength to a finished coating 
due to its own high hardness (7,5 - 8.0 in the Mohs hardness 
scale). Furthermore, the capacity of zinc spinel to form solid 
solutions and its wide isomorphism can be used for the syn¬ 
thesis of colored coatings in the presence of chromium oxide, 
and color centers based on ZnO • A1 2 0 3 will be resistant to 
the effect of aggressive glaze melts. 

Thus, by predicting the reaction products in the system 
considered at the temperatures of the melt formation, the fol¬ 
lowing has been established: 

- all phase transformations in composites of series 3 and 
Z occur at temperatures below 1500 K, which makes it possi¬ 
ble to perform heat treatment of coatings at a temperature not 
higher than the specified one; 

- analyzing the final phase composition of heat-treated 
products, it is necessary to take into account the firing dura¬ 
tion of non-fritted glazes (36 - 48 h), due to which such 
compounds as helenite, sillimanite, and forsterite will most 
probably pass into the melt, and others, such as gahnite, 
mullite, and diopside, on the contrary, will be crystallized 
from the melt; 

- with respect to obtaining colored coatings, the most 
suitable interval is 1400 - 1500 K, since magnesium or zinc 
spinel most probably are formed within this interval, and fa¬ 
vorable conditions are developed for the formation of solid 
solutions on their basis; 
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- a color-bearing phase in the optimum coatings of series 3 
will presumably be a solid solution of iron, MgO • Al-,0 3 , and 
BaO • Fe 2 0 3 , although it is hard to predict the color of coating; 

- the optimum coatings of series Z most probably will 
get their tint from spinelide of complex composition 
Zn(Al, Cr),0 4 , and the glaze of composition Z 21 is likely to 
have the most intense color. Synthesized compositions of 
non-fritted glazes of the series considered fully corroborated 
the optimum heat treatment interval found in the calculation, 
namely, 1503 -1553 K. High-quality tinted (mustard and 
inol color) coatings with good luster, thermal resistance, and 
microhardness were developed due to the synthesis of 
spinelide compounds. 

The phase composition of coatings was determined using 
x-ray phase analysis (Figs. 3 and 4), which corroborates the 
calculation validity with respect to both quantitative and 
qualitative phase content. 

Thus, the use of the Karapet’yants method together with 
the Hooke - Jeeves method and the developed method for 
setting initial approximations to calculate the system of non¬ 
linear equations has made is possible to make a modeling 


calculation of equilibrium compositions obtained at various 
temperatures from various batch mixtures. 

A possible mechanism of phase formation for varying 
temperatures of the system has been identified and the ap¬ 
proximate content of reaction products and ensuing proper¬ 
ties of the finished coating have been estimated. An optimal 
heat treatment procedure has been selected to obtain opti¬ 
mum service parameters in coated products. 
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